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SECTION  I 
INTRODUCTION 


This  report  describes  research  which  was  directed  at  the 
development  of  an  orthotropic  plate  finite  element  for  the 
analysis  of  plate  and  shell-like  structures  which  exhibit 
coupling  between  extension  and  bending.  The  element  is 
especially  useful  in  the  analysis  of  structures  which  are 
fabricated  from  laminated  composite  materials.  The  report  is 
written  so  that  it  would  describe,  for  the  finite  element 
expert,  the  analytical  techniques  utilized  in  the  development 
of  the  element  ana  also  would  be  of  use  to  a  program  "user" 
not  having  the  overall  expertise  of  an  expert. 

The  notation  and  the  methods  used  for  the  definition  of 
element  material  properties  have  been  chosen  as  a  result  of 
a  careful  survey  of  the  literature  on  composite  materials. 

The  notation  and  definitions  chosen  are  considered  to  be  in¬ 
dustry  standard  and  are  best  summarized  in  reference  3,  which 
is  rapidly  becoming  a  standard  text  for  the  analysis  of  com¬ 
posite  materials. 

An  industry  standard  computer  program  SAP  IV  ^ ^  was 
selected  as  host  program  to  accept  the  new  composite  plate 
finite  element.  The  SAP  IV  Finite  Element  Computer  Program 
was  designed  to  easily  accept  new  elements  into  its  element 
library.  The  new  element  muot  be  self  contained  since  the 
general  philosophy  and  program  structure  is  "overlayed"  into 
the  computer.  The  laminate  composite  plate  element  is  the  new 
element  to  be  integrated  into  the  element  library.  Called 
TYPE  9,  the  new  element  is  similar  to  the  SAP  element  TYPE  6 
in  both  description  and  input.  The  main  difference  is  that 
element  TYPE  9  has  the  ability  to  describe  the  effects  of 
coupling  between  in-plane  extension  and  out-of-plane  bending. 
Element  TYPE  9  is  a  quadrilateral  element  and  is  formulated 
from  quadrilateral  shape  functions  rather  than  from  four 
triangles  as  in  TYPE  6.  Also,  TYPE  9  allows  material  directions 


1 


o  be  arbitrary  for  ease  of  material  input  descriptions.  The 
lament  is  modelled  after  the  structure  of  element  TYPE  6; 
herefore  element  TYPE  9  can  degenerate  to  element  TYPE  6. 
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SEC  T  I  ON  1 1 

LAMINATE  COMPOSITE  FuAT  PLATE  ELEMENT 


The  laminate  composite  flax  plate  element  is  based  on  thin 
plate  theory  with  the  exclusion  of  transverse  their  deforma¬ 
tions.  The  following  sections  describe  the  basic  f jrmulat ion 
of  the  finite  element. 

1.  ELEMENT  POTENTIAL  ENERGY  FUNCTIONAL 

The  principle  of  minimum  potential  energy  furnishes  a 
variational  basis  for  the  direct  rormulation  of  the  element 
stiffness  equations  and  loading  functions.  The  potential 
energy  of  the  element  is  formed  from  the  sum  of  strain  energy 
(U^)  and  the  potential  of  all  applied  loads  (V  );  i.e., 


The  principle  can  be  stated  as  follows:  Among  all  the  displace¬ 
ment  functions  of  admissable  form*  those  that  satisfy  the 
element  equilibrium  conditions  make  the  potential  energy  func¬ 
tional  obtain  a  stationary  value.  Thus, 


6rr  =  6U  +  5V  -  o  ( f ) 

P  p  ? 

where  5  is  the  first  variational  operator. 

It  can  be  shown  that 

6U=/ai6edV  ( 3 ) 

P  >  - 

where  o  is  a  vector  of  stress  components, 

t  the  corresponding  vector  of  strain  components  and 

V  the  volume  of  the  element. 

All  vectors  will  be  underscored  with  a  straight  bar  and 

matrices  will  be  underscored  with  a  tilda.  The  super- 
T 

script  of  the  vectors  and  matrices  designates  the 
matrix  is  transposed. 


Note : 


The  corresponding  first  variation  of  the  potential  forces 
become s 

5Vp  =  -  f  bT  5u  dV  -  f  £T  6u  dS  (4) 

J  v  ^  s 


where  b 

repre 

seats  the 

element 

body  forces  , 

x  L>  vl 

vector  or 

surface 

tractions  applied  on 

3ur: ace  S,  and 

u 

is  a 

vector  of 

element 

displacements . 

here  that  the  surface  traction  integral  can  be  used  to  include 
the  point  concentrated  forces  on  the  boundaries  of  the  element. 

The  elements  of  the  strain  potentials  of  equation  (3)  will 
eventually  lead  to  the  element  stiffness  and  initial  load  vec¬ 
tors  and  the  elements  of  the  applied  load  potential  of  equation 
(O  will  produce  the  various  element  vectors. 

2.  QUADRILATERAL  SHAPE  FUNCTIONS 

The  element  formulation  is  a  geometrically  linear  quadri¬ 
lateral  containing  the  four  corner  nodes  as  shown  in  FIGURE  1. 


FIGURE  1  -  Quadrilateral  Element  Geometry 
a.  Geometric  Shape  Functions 

The  element  shown  in  FIGURE  1  is  described  in  the  local 
coordinate  of  the  element  and  all  material  reference  is  made 
with  respect  to  the  element  local  x  axis.  The  element  area 
domain  can  be  described  by  using  a  polynomial  as 


4 


(5) 


T 

x  =  £  6 

T 

y  =  i  6 

where  x  and  y  are  the  local  coordinates  as  in  the  element 
domain , 


r  Li,  r,  s,  rsj  (6; 

the  row  vector  of  polynomial  coefficients  jS,  a  vector  of 
generalized  coefficients,  and  r,s,  the  element  natural 
coordinates . 

The  generalized  coefficients  can  be  solved  for  by  evaluating 
the  polynomials  at  the  vertices  of  the  element.  Therefore, 

x  =  HT  x  (7) 

T 

y  =  h  £ 

T 

where  H  contains  the  terms  of  the  shape  function  of  the 
element;  x  and  y  are  vectors  containing  the  element  ver¬ 
tices  as , 


x  ^  =  Lx-^ ,  x2,  x  ^  5  x4J  (8) 

T 

y  =  Lyx,  y2,  y3,  y4J 

The  terms  of  the  shape  function  can  be  described  by 

hi  =  —  (1  +  r^r)  (1  +  s.s)  (9) 


where 


£iT  =  L-l ,  1,  i,  -1J  (1C) 

siT  =  L-l,  -1,  1,  1J 

define  the  natural  coordinates  ot  the  elements. 


The  mapping  of  the 
can  be  obtained  by 


element  geometry  and  displacement  functions 
defining  the  Jacobian  transformation  as 


3  -  J  ^ 

-n  -  -£ 


(11) 


where  J  is  the  Jacobian  matrix  defined  as 


and 


(12) 


(13) 


are  the  first  derivative  operators  in  the  natural  and  local 
reference  frames  respectively. 

Note:  The  " , ,r  subscript  implies  "partial  differentiation  with 
respect  to".  The  inverse  transformation  is  obtained  by 


3 


Z 


(15) 


where  J*  is  the  determinant  of  the  Jacobian  matrix  given  as 


J* 


’r 


(16) 


and , 


6 


It  will  be  necessary  to  obtain  second  derivates  in  the 
local  reference;  therefore 


4 


e  a  +  f  a: 

-  — n  -  — n 


(1 


represents  the  second  partial  operator  in  the  local  reference 
given  as 


8x9y 


and  the  natural  set  as 


32 


The  E  matrix  is  defined  as 


where 


e  .  . 


£  2  gT 

J*  §i 


-n 


(22) 


with  gt  being  the  ith  row  partition  out  of  the  G  matrix  and 


„  *  1  < 
3  =  3  -  4*  J 

-n  -n  J  *  -n 


(23) 


The  F  matrix  is  defined  as 


where 


with 


T  , 

hi 

f  t 
-22 

m 

f  x 

-12  J 


(24) 
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J*2  Si  li 


(25) 


gi 


gil 

gi2  ° 

o 

gil  gi2 

(26) 


the  elements  being  of  the  G  matrix, 
b.  In-Plane  Displacement  Shape  Function 

The  plate  element  is  assumed  to  have  in-plane  deformations; 
therefore  the  variation  of  x  and  y  displacements,  u  and  v  res¬ 
pectively,  can  be  expressed  using  the  same  shape  functions  for 
the  geometry  as  in  the  previous  section.  Then 


u 


V 


(27) 


8 


the  element 


are  the  domain  displacements  of 


where 

and 


T  . 

L u,  u„ 

u„  u  J 

u 

1  2 

o 

4 

m 

i 

< 

> 

> 
— i 

Vo 

v,  J 

V 

1  L 

3 

4 

(28) 


contain  the  in-plane  model  displacements.  Therefore,  it  is 
assumed  that  in-plane  displacements  vary  linearly  within  the 
element . 

c.  Transverse  Displacement  Shape  Function 
The  plate  element  defined  by  thin  plate  theory  must  have 
a  transverse  shape  function  to  allow  for  proper  bending. 
Therefore  it  is  assumed  that  the  shape  function  polynomial  is 

J  ,  ,  2  232  233  3, 

$  =  Ll,r,s,r  ,rs,s  ,r  ,r  s,rs  ,s  ,r  s,rs  J 

(29) 

The  natural  degrees  of  freedom  allowed  per  node  for  bending  are 

2oi  =  Lw  w’r  w’sJ  i  (30) 


Letting 


q0  =  f  §  (31) 

where  g-  is  the  full  set  of  degrees  of  freedom  in  the 

natural  reference  of  the  element,  i p  is  the  matrix 
$  evaluated  at  the  nodes  and  defined  in  reference 
2  and  6  is  the  set  of  generalized  nodal  coefficients. 
Then  the  transverse  displacement  w  becomes 

w  =  $n  *  ^0  -  90  (j2) 

where  ij>  ^  is  the  inverse  of  ip  and 

il  i:,  i  lie  I  r*. in ■ ; vi* r1: i* *  uh.ipo  function,  both  o!  which 
1 1  rt  •  Ji  i  i  1 1*  *<  1  .  n  r< '  I  <  •  r ’  1 1<  ■<  ‘  V  • 
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The  transverse  displacement  shape  functions  are  defined  in  the 
natural  reference  of  the  element  to  allow  for  ease  of  develop¬ 
ment  since  second  derivatives  must  be  taken.  Since  a  global 
transformation  by  node  is  to  be  performed  at  a  later  stage, 
the  transformation  by  node  from  the  natural  to  local  coordi¬ 
nate  can  be  made. 

3.  PLATE  STRAIN  FUNCTIONS 

The  classical  assumptions  of  linear  thin  plate  theory  are 
made,  essentially  reducing  the  three-dimensional  equations  of 
elacticity  to  a  two-dimensional  set  of  plane  stress  equations. 
For  the  elastic  continuum  of  the  plate,  the  following  assump¬ 
tions  are  made: 

•The  thickness  (h)  is  small  compared  to  the  dimensions  of 
the  plate  in  the  x  and  y  directions. 

•A  line  element  through  the  thickness  remains  normal  to  the 
mid-plane  surface  under  all  states  of  deformation,  inde¬ 
pendent  of  its  translation  or  rotation. 

•The  plate  can  be  isotropic,  orthotropic  or  comprised  of  a 
number  of  orthotropic  laminae,  where  each  lamina  obeys 
Hooke’s  law. 

•The  displacements  u,  v,  and  w  in  the  x,  y,  and  z  direc¬ 
tions  respectively,  are  small  when  compared  to  the  plate 
thickness . 

•The  reference  axis  is  taken  as  the  middle  of  the  plate  at 
h/2,  h  being  'the  total  plate  thickness. 

•The  normal  strain  in  the  z-direction  is  assumed  to  be 
zero,  giving 


£ 


Z 


W,. 


0; 


therefore,  the  lateral  deflection  is  given  by. 


w  -  w(x,y)  . 
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•  St.  Variant1  s  principle  applies.  That  is,  local  deformation 
occurs  in  the  area  of  applied  loads  while  at  distances 
away  from  the  load,  the  deformation  state  is  not  grossly 
affected . 

•Transverse  shear  deformations  are  neglected; 

Y  =  Y  =  0  . 
x  z  yz 

•Displacements  are  linear  such  that 

u  -  Uq  (x ,  y )  -  zw,x  and  v  =  vQ(.x,y)  -  zw,y, 

where  w,  =  -6  .  w,  -  9  ,  and  ur  and  are  the  in-plane 
5  x  y  y  x  o  o 

displacements  of  the  middle  surface.  The  rotations  about 

the  x  and  y  axes  are  given  by  ©  and  Q  ,  respectively. 

x  y 

a.  Midplane  Strain  and  Curvatures 

The  mechanical  strains  associated  with  plate  stretching 
and  bending  can  be  written  as 

£re+z<  (33) 

-o 

where  the  mid-plane  strains  are 

(  3^  ) 


the  plate  curvatures  are 

xx  1 

(35) 

yy  ( 

'  xy 


±  L 


with  u  and  v  being  the  in-plane  displacements  and  w,  the  trans 
verse  displacements.  The  thermal  strains  can  be  written  as 


e.  =  a  T  +  z  a  T 

-t  -  o  -  g 


(36 


where  a  is 
relative  to 
erature  dif 
through  the 


a  vector  of  thermal  expansion  coefficients 
mid-plane  strains,  Tq  the  element  mean  temp- 
ierence  and  the  element  thermal  gradient 
plate  thickness. 


b.  Strain  Displacement  Functions 

The  connection  between  strain  and  displacement  is  made 
realizing  that  the  in-plane  and  transverse  displacements  have 
been  made  relative  to  a  set  of  nodal  displacements.  Equation 
(33)  can  be  written  as 


£  =  ?I  Si  +  z?0  20 


(37 


where  is  the  in-plane  strain-displacements  relative  to 
(the  in-plane  nodal  displacements); 

is  the  transverse  strain-displacements  relative  to 
(the  natural  transverse  nodal  displacements). 

The  in-plane  aisplacements  are 


Si 


(38 


and  the  3.  becomes 


where 


..T 

n 

~n 


(40) 


The  out  of  plane  displacements  are 


So 


3o2 


^03 


(41) 


i3°4J 


where  the  sub-elements  of  the  partition  are  defined  by 
equation  (30)  and  becomes 


i  [e  R;  +  f  flj  ] 

~3  ~  ~n  ~  -nn 


(h2) 


where 


ZZ 


[-1  0 

0  -1 
0  0 


(43) 


5  FT 

-n  - 


(  4  4  ) 


ana 


-T 

H 

-nn 


2  — ' T 

2„  5 


(4b) 


Note:  In  equation  (37)  the  z  variable,  which  is  the  plate 1 s 

normal  coordinate,  is  maintained  distinctly  since  it  is 
independent  of  the  in-plane  variables.  Later,  when  the 
strain  energy  is  formed,  the  z  variable  will  integrate 
through  the  thickness  and  merge  into  material  property 
matrices . 
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PLANE  STRESS  COMPONENTS 


The  stress  components  for  a  thin  plate  can  be  written  in 
vector  form  as 

o [  *  C  (e  -  eT>  (46) 

where  C  is  the  material  matrix  described  in  the  plate  local 
axes  and  is  expressed  as 

£  =  RT  C  R  (47) 

^  C  *  ~  E 

with  C  being  the  material  matrix  in  the  principal  material 

m 

directions  of  the  fibers  and  R1  being  the  strain  trans¬ 
formation  matrix  from  element  focal  coordinates  to  prin¬ 
cipal  fiber  directions.  The  elements  of  equation  (47) 
are  found  in  reference  2. 

The  elements  of  the  thermal  strain  involving  the  thermal  coef¬ 
ficients  are  defined  as 


ot  =  Rr  a  . 


(48) 


Once  the  material  matrix  is  defined,  the  elements  of  the  material 
...Lasticity  matrix  can  be  defined  as  in  the  following  section. 

i.  MATERIAL  ELASTICITY  MATRICES 

The  material  coefficients  are  defined  with  the  use  of 
equation  (3)  in  a  slightly  different  form: 


e  dV 


(49) 


Since  the  local  z  dimension  is  small  compared  to 
plate  dimensions,  it  is  convenient  to  define  the 
tants  and  moment  resultants  as 


the  x  and  v 
stress  resul- 


=  ft  a- 


d  z 


(SO) 


4 


^na 


'hen 


*  ;  L 


a  z  d  z 


a  new  stress-strain  matrix  can  be  defined  as 


(Si) 


where 


o  - 


M 


'A  31  l  e  j 

-  -G 


|_bT  d  |  ts 


(S_ 


Letting 


and 


1  fj 

*  /  s 


C  dz 


C  zdz 


b1 


C  z  dz 


A^,  Bip 

B_  Dm 
-T  -T 


(53) 

(54 ) 

(Si ) 


(56) 


(57) 


(58; 


equation  (49)  can  be  written  as 


X  2* £  i 


~T  p  — 

£  r  £  dA 

m  aA 


/ 1 


T  E  J'oi 


-T,TjdA 
(59) 
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where 


A  T 

■  A 

c 

CL 

dz 

(60) 

-T 

-A 

c 

ot 

zdz 

2 

(61) 

-T 

"  A 

c 

a 

z  dz  . 

(62) 

The  material  A,  B,  D  matrices  and  the  thermal  load  coef¬ 
ficients  A^,  B  and  D  can  be  related  to  laminar  material  by  posi- 
tion  t  in  the  material  build-up  as 


L 

A  = 

C.  (t.-t.  , ) 

(63) 

i  =  l 

~  X  -L  X  ^  X 

B  = 

1/2 

t 

(64) 

D  = 

1/3 

Ix  Si  <t l-  tf.,) 

(65) 

L 

A„  = 

—  i 

l 

1  =  1 

Si  Si  (ti  -  ‘i-l’ 

(66) 

-T 

1/2 

T  c.  o.  (t?  -  t?  .) 

i^l  1  1"1 

(67) 

-T 

1/3 

t  Si  Si.  <*!  - 
1=1 

(68) 

where  the  subscript  "i"  implies  coefficient  evaluation  at 
laminae  level  Mi"  and  L  is  the  total  number  of  fiber  lamina 
levels . 

6.  COORDINATE  TRANSFORMATIONS 

The  element  information  is  initially  determined  in  the 
natural  coordinates  of  the  plate  since  it  is  quite  easy  to  expres 
all  loading  and  stiffness  information  in  that  reference.  Ulti¬ 
mately  the  information  must  be  transformed  to  local  coordinates 
(x,y,z)  and  also  to  global  coordinates  (X,Y,Z).  The  following 
sections  describe  the  transformations. 


u.  Natural  to  Lo^ai  Trur.u: ormat ion 

The  natural  coordinate  variables  per  node  are  defined  a 

/  \ 


J*  \ 

'•*ni  v  w  ’  r  /  i 

!  w’s 

where  6 is  the  rotational  degree  of  freedom  normal  to 
plate  at  node  "i". 

The  transformation  matrix  required  becomes 


where 


luvwu  G  GJ 
x  y  z 


with 


lw'r)  j°x 

f- '  "  Si  k 


Noting  that  the  rotation  degrees  of  freedom  are  defined 


.n«=  v  matrix  cccomes 


The  complete  natural  to  local  transformation  therefore  becomes 


~  nj£l 


~n£  2 


b.  Local  to  Global  Transformations 

The  local  to  global  transformation  quantities  are  somewhat 
more  difficult  to  obtain  since  the  transformation  involves  the 
local  coordinates  of  a  quadrilateral  element.  Obviously  only  3 
points  define  a  plane;  the  fourth  point  of  the  quadrilateral  is 
unnecessary.  The  fourth  point,  however,  may  not  lie  in  the  same 
plane  as  the  other  three  points.  Therefore  local  transformations 
by  node  are  determined  and  are  averaged  to  obtain  a  general 
transformation  used  in  determining  the  element  coordinates  and 
m  transforming  element  matrices  where  applicable. 

Defining  the  nodes  of  Figure  1  as  i,  j,  k  and  1  and  allowing 
this  sequence  to  permute,  the  element  normal  coordinate  at  node 
"f  which  also  permutes,  is 

V  .  .  X  V  . 

S  _ Tht  (76) 

^  zi  IV.  .  X  V.  J  w  ; 

-:i  -ii 

where  the  "X"  symbols  denote  a  cross  product  of  two  vectors, 

and  V . ■  implies 
-li  F 
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(77) 


with  X,  Y  and  Z  being  global  coordinates. 

A  global  material  reference  is  defined  as  a  global  vector 
defined  as  input  which  locates  the  general  local  "x"  of  all 
elements  referred  to  that  vector.  All  material  properties  are 
defined  relative  to  this  x  coordinate  which  becomes  the  ele¬ 
ment's  local  x  coordinate.  The  element's  local  y  coordinate 
at  each  node  can  be  calculated  as 


e  .  X  V 
21  -m 

ITT-] 


zi 


(70 


Finally,  the  local  x  coordinate  at  each  node  is  determined  as 

(79) 


§  .  X  g  . 
yi  zi 


xi 


!e  .  X  §1^ 
t  yi  zi 1 


The  local  to  global  transformation  at  node  "i"  becomes 


L  • 
-  1 


reTl 

X  I 


1 

ey 


(60 


where  the  direction  cosines  of  each  coordinate  are  placed 
in  row  order  in  the  transformation  matrix. 

The  average  transformation  used  to  determine  the  local 
coordinates  is  obtained  by  first  averaging  the  nodal  normals  as 


e 

z 


(81) 


and  substituting  into  equations  (78,79  and  80)  to  produce  an 
average  t.  . 


19 


The  global  degrees  of  freedom  can  be  defined  by  node  as 


'u  s 

V 

w 

6, 

e. 


(82) 


v  / 

where  U,  V  and  W  correspond  to  global  displacements  rela¬ 
tive  to  X,  Y  and  Z  respectively  and  0fs  corresponds  to 
global  rotations  about  X,  Y  and  Z  respectively. 

The  complete  transformation  becomes 


where 


~£gi 


■lg2 


£g4 


(83) 


:zgi 


(84) 


and  and  are  t^e  complete  list  of  degrees  of  freedom 
per  element. 

If  the  quadrilateral  element  is  perfectly  flat  in  its  space, 
then  the  •  becomes  exactly  tn  .  If  it  is  not,  then  the 
element  space  appears  as  a  curved  space.  This  effect  should 
allow  the  element  to  behave  as  a  shallow  shell. 


c.  Natural  to  Global  Transformat  ion 

The  element  stiffness  will  be  transformed  from  natural  to 
global  coordinates  directly.  Therefore  that  transformat  ion 
becomes 


20 


r 


i 


£ 


n 


(85) 


where 


-ng 


(86) 


Since  both  original  transformation  matrices  are  partitioned 
diagonally,  the  nodal  transformation  matrix  is 


t  .  =  t  .  t  . 

~ngi  ~n£i  ^£gi 


(87) 


which  contains  many  off-diagonal  zeroes.  Therefore  it  is  con¬ 
venient  to  define  transformations  by  displacement  and  rotation 
degrees  of  freedom.  That  is,  the  displacement  transformation 
at  node  i  is 


:Di 


(88) 


and  the  rotation  transformation  at  node  i  is 


(8S) 


This  modification  will  save  transformation  operations  later  on. 


7.  ELEMENT  STIFFNESS  MATRIX 

The  element  stiffness  is  easily  defined  in  the  natural  coor 
dinates  of  the  plate,  given  the  strain  displacement  functions. 
Once  this  stiffness  is  determined,  it  can  be  augmented  with  a 
scaffolding  or  artificial  torsional  stiffness  for  the  plate’s 
normal  degrees  of  freedom.  Finally,  this  stiffness  can  be 
transformed  to  local  and  then  to  global  coordinates  for  assembly 
into  a  master  stiffness  matrix.  The  following  details  the  above 
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a.  Plate  Element  Stiffness  in  Natural  Coordinates 
Defining  a  new  set  of  degrees  of  freedom  as 

5-  •{£!• 

the  strain  components,  as  defined  by  equation  (58)  become 


B 


~n 


^n 


(91) 


where 


(92) 


The  second  integral  of  equation  (59)  can  be  used  to  define  the 
element  stiffness  in  local  coordinates  as 


Em  e 
~m  - 


dA 


qT  K  q 
-n  ~n  ~n 


where 


R 


f  BT  E  B  dA  . 
J  ~n  ~m  ~n 

A 


(93) 


(94) 


For  convenience  of  computation,  the  material  matrix  Em  in  equa¬ 
tion  (94)  is  Cholesky  factored  as 


T 

E  =  U  U 
~m  ~ 


(95) 


where  U  is  a  upper  tri-diagonal  factoring  matrix. 

This  allows  equation  (94)  to  be  written  in  a  more  efficient 
form  as 

K  -  f  (U  B  )T  (U  B  )  dA  (96) 

A 

which  allows  the  triple  matrix  product  to  be  replaced  by  a  sim¬ 
pler  transpose  symmetric  product.  This  process  is  especially 


2  2 


efficient  since  the  Choi u sky  factoring  is  performed  (at  most) 
one e  per  element.  Also,  numerical  integration  is  to  be  per¬ 
formed.  Savings  will  occur  at  each  integration  point  after 
the  first. 

b.  Artificial  Torsional  Stiffness 

The  flat  plate  theory  does  not  have  any  mechanism  to  di¬ 
rectly  include  twisting  of  the  plate  normal  to  the  plate  sur¬ 
face.  Therefore,  if  two  coplanar  elements  are  assembled  at  a 
common  node,  a  singular  stiffness  exists.  To  avoid  this,  an 
artificial  or  scaffolding  stiffness  is  added  to  the  normal 
rotational  degree  of  freedom  G  .  There  is  no  change  in  the 
system  equilibrium.  For  convenience,  this  is  performed  at  all 
nodes  of  the  element  since  it  would  be  difficult  to  determine 
coplanar  effects  in  general.  Tnis  does  change  the  overall 
element  equilibrium.  If  the  amount  of  artificial  stiffness  is 
Kept  small  and  the  local  rotational  stiffness  effects  are  in 
equilibrium,  then  the  error  can  be  minimized.  Defining  a  vec¬ 
tor  of  normal  rotations  at  the  nodes  as 

(O 


the  artificial  torsional  stiffness  matrix  relative 


:o  6  becomes 


-1 

-1 


where  f  is  an  input  scaling  factor  which  can  vary  as 


0  <  f  <  1 


(99) 


and  can  be  set  in  l.E-8  increments,  C  is  an  artificial 
coefficient  estimated  from  element  bending  stiffness 
coefficients  and  element  area;  i.e., 

C  =  MIN  (D( 1 , 1 ) ,  D ( 2 , 2 ) )  *  AREA  (100) 

with  the  D's  defined  in  equation  (55). 
c.  Natural  Stiffness  Matrix 

The  degrees  of  freedom  q  and  6  defined  by  equations  (90) 

—  n  z 

and  (97)  can  be  merged  to  degrees  of  freedom  gn  described  in 
equation  (69).  This  requires  the  re-ordering  of  stiffness  coef¬ 
ficients  of  equations  (96)  and  (98)  to  produce  a  natural  stiff¬ 
ness  matrix: 


K 


~n 


MERGING 

REORDERING 


(101) 


relative  to  q  . 

-*n 

d.  Global  Stiffness  Matrix 

The  global  stiffness  matrix  K  is  formed  by  transforming 
K  from  natural  coordinates  to  global  using  equation  (85).  The 
transformation  is  formed  using  equation  (1)  realizing  that  the 
strain  energy  U  is  invariant  relative  to  any  coordinate 
reference.  Therefore, 


1  T 

U  =  4  q  K 
p  2  -n  ^n 


Using  equation  (85)  produces 

T  ' 

_  m  -L 


2  9g 


K  q 
~g  -g 


(102  ) 


K 


<tg 


K 


(103) 


Tin  triple  rn.  1 1  r  •  i  j  r  «lui  I  i  mp  I  i  i  •<  I  in  <  <pi.i  I.  i  <  hi  (  1  0  O  i :  q  u  i  tr¬ 
ine  f  f  ic  lent ,  e  u  pec  i  .1 1 1  y  l;  i  nee  Y  i..-.  highly  <1  i agonal  .  Efficiency 
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sub-.natr  ices  labelled  X 
to  8 .  Then 


uy  part  it  toning 

n  .  g 
.  ana  n -  .  wnere 
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K  and  K  into 
i  and  i  range 


3X3 
rom  1 


(104) 


wnere  t. 
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9  > 

4, 


5, 


matrix  relates  to  equation  (38)  when  i  equal 
7  ana  relates  to  equation  (83)  when  i  equals 


Additional  eri'ieioncy  is  obtained  wnen  ttc  triple  matrix  product 

is  performed  such  that  the  right  portion  matrix  multiply  is 

first  formed  ana  positioned  back  into  K.1?.  Then,  the  left 

-i  j 

multiply  is  formed  with  the  resulting  product  ana  placed  into 

rr 

X  .  V  . 

~  -  J 


c .  NUMERICAL  INTEGRATION  OF  AREA  FUNCTIONS 

The  elements  in  the  matrix  of  equations  (94)  and  (96)  are 
very  difficult  to  integrate  exactly,  therefore  approximate 
numerical  integration  can  be  performed  with  sufficient  accuracy 
for  justification.  Gauss-Legendre  Numerical  Quadrature  has 
oeen  selected  to  perform  the  integration  of  the  stiffness  coef¬ 
ficients  as  well  as  other  area  functions.  Equation  (64)  can  be 
rewritten  and  transformed  relative  to  variables,  and  limits  of 


prat: 


A 

-  n 


^T 

B 

~  n 


(r 


^  ^  tjnb'n 


( r  ,  s  ) 


eras  . 


(  _  6  ) 


An  element  of  chis  matrix  can  be  written  as 

i:  r 


K  .  . 
i  : 


-1 


:  .  ( r ,  s  )  drds 


Clue) 


where 


:  ;  ,  s  >  -  E  bJ  (k,i)  E  i:rW,{)  :-rU,i) 

■  '  -n  P:  n 


(-07) 
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The  stiftness  coefficient  can  then  be  approximated  as 

%  fij  <w  «k  w<  a 

where  n„ ,  n0  are  the  number  of  Legrendre  root  evaluation 

X  4 

points  in  the  r,s  directions  respectively, 
r.^  ana  s0  are  the  roots  of  the  Legrendre  polynomial, 
and  are  the  appropriate  Gauss  weighting  factors. 


For 


r  or 


£  £ 

rk  =  sk  MT  ’VT 
Wk  =1,1 


(109) 


n  =  3 , 


yr  VF 

5  ,  0,  5 

5  8  5 

9  ,  9,  9 


(110) 


3.  ELEMENT  MASS  MATRIX 

A  consistent  mass  matrix  relative  to  the  in-plane  variables 
in  one  coordinate  can  be  written  as 


L 


p  t  H  H  dA 


(111) 


and  a  corresponding  lumped  mass  matrix  can  be  formed  by  summing 
the  rows  of  the  consistent  mass  matrix  as 


M 


l 


p  t  f  HT  dA 


assuming  t  and  p  constant  and  realizing 


(112) 


z.  Hi 


i. 


26 


he  co^cner.tL  n;  are  applied  to  the  translatory  degree^ 
reedem  in  ail  ^ODal  dire  cions,  per  node.  A  rotary  inerti 


ciiect  can  be  included  by  an  approximation; 


:nt-  gioDai  iumpcu  mass  vector  can  oe  formed: 


m2  i 


: 


where  i  is  a  o  x  i  unit  vector. 

1C.  ELEMENT  LOAD  VECTORS 

The  element  load  vectors  are  established  from  element  p 
erties  such  as  material  constants,  temperatures,  pressure,  m 
area  and  acceleration  constants.  The  following  sections  ces 
cribe  the  load  vectors  •developed. 

a.  Thermal  Lvid  Vector 

The  there  integral  of  equation  (59)  can  be  usee  tv  def. 
the  clement  t normal  load  vector  using  equation  (SI): 


7  T  dA  =  a; 


B‘  dA  i:,„  > 

-1  (xg) 


mere  r  ore , 


B  dA 
-  n 


is  trie  thermal  vector  relative  to  the  natural  coord  i  nu  tot 


.)  i  :  I 


i  i  Vi  *• '  I  ■  >r  ■  !  *  » »  ,iii<  I  In  ij:,  : .  t  r  ■.  i  i 


r 


a,r 


p  J  O' 

?t  m 
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rhe  global  thermal  vector  can  be  determineu  by  applying  the 
natural  to  global  transformation. 


rP  rp 

T  "  F x 
-ng  -n 


(116  ) 


b.  Pressure  Load  Vector 

The  pressure  load,  by  node,  is  formed  for  tho  shape  terms 
associated  only  with  the  geometry  and  is  applies  in  the  2  (plate 
normal)  direction.  Therefore, 


f  p  H  dA 

J  A 


(119) 


A  pressure  load  vector  can  be  formed  relative  to  local  cocrdi- 


na t e s  -  as 


/  \ 
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Ei!  < 
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(120) 


?he  global  pressure  load  can  be  formed  as 


7  1  pP ' 
i£g  i£ 


(121 ) 


which  transforms  the  normal  traction  into  a  global  traction, 
c.  Constant  Acceleration  Load  Vector 

The  acceleration  vector  can  be  computed  from  the  mam,  wet 02 
defined  in  equation  (114);  i.e.. 


a  m 


-g 


(  L  2  ?  ) 
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MODIFICATIONS  TO  SAP  IV  COMPUTER  PROGRAM 

The  SAP  TV  computer  program  wcii  modified  to  accept  the  new 
opposite  element.  The  following  sections  describe  changes  to 
he  existing  progiMin  as  well  as  new  routines  of  the  composite 
Ic.  tt-  e  Leniun  t  . 

SAP  IV  STRUCTURE 

The  main  changes  to  the  SAP  IV  program  occur  in  the  elemen 
ibrary  control.  They  occur  in  the  clement  generation  portion 
f  the  program.  Figures  2  and  3  depict  the  major  routines  used 
y  SAP  to  perform  a  static  analysis  of  its  elements.  The  FI.TYP 
satire  calls  a  new  routine  called  PLATE  *  which  is  usee  to  cd.  i 
VERLAY  (8).  (stress  recovery)  Once  OVERLAY  (8)  is  called,  a 
all  is  made  to  CP LATE,  the  main  routine  of  the  composite  plate 
i.ilte  element.  Figures  4,  5,  6,  and  7  describe  the  flow  of 
-on tines  used  by  CPLATE. 


. .  w 


DESCRIPTION 


Input  title  ana  master  control  card 
Input  nodes 

Write  8  ID,  X,  Y,  Z,  T 

Call  element  type  (new  composite 
plate  element  called  here) 

WRITE  1  Element  stress  matrices 
for  recovery 

WRITE  2  Element  stiffness  matrices , 
mass  and  load  vectors 

Input  Concentrated  I oads/masses 
WRITE  3  Concentrated  loads/ 
masses  for  assembly 


READ  2  Read  stiffness  ror  asoembiv 
WRITE  4  Assembled  Stiffness  and 
Loads  in  blocks 
WRITE  3  Assembled  Mass  Vector 


Call  appropriate  solution  type 
0  Static 

1  Nodal  Extraction 

2  Force  Response  Analysis 

3  Response  Spectrum  Analysis 

4  Direct  Integration 


FIGURE  2  BASIC  PROGRAM  FLOW 


Static  solution  option 

Solution  or  displacement  equation 
READ  4  Read  Stiffness  and  ioado 
for  decomposition 
WRIT;!  2  Solve  for  displacements 

READ  8  Read  ID  i dent  if icat ion 
matrix 

READ  2  Recover  displacements 
Print  displacements 

READ  1  Element  stress  matrices 
Stress  recovery  and  Printing 


SOLUTION  AND  RECOVERY  STAGE 


DESCRIPTION 


Composite  Plate  routine 

input  material  Properties  in 
tables 

Input  global  material  vectors 


Input  load  factors  TLO  as  a 
^unction  of  thermal,  pressure 
and  constant  acceleration 


Input  element  definition 


Form  ABD  matrices  and  calcu¬ 
late  ov  thermal  stress  vector. 
ModifyAE  =  U*'  U 


Determine  local  to  global  trans¬ 
formation  and  locai  element 
coordinates . 

OROSN  -  forms  vector  cross  product 
STRANS  -  special  transpose 
KG  I.  -  used  to  transform  global  to 
local  coordinates 


bet  element  K  , 
.  -n 

to  zero. 


M.^  matrices 


.  ’  I-fo  h 


4  FLOW  DIAGRAM  TOR 
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v  ailm  a  description  o.  the  routines  used 
_ it  r 1 a t e  overlay. 

•  i  ue  •  ■  *1  lea  by  CFLATI!  to  fill  tne  A,  h  arc  D  02 
S3,  :  and  55  material  matrices  :  rom  the  input 
•  ient 0  arrays.  equations  CO,  Cl,  and  62  are 

-  :  it*--  tin  thermal  load  vectors,  and  thermal 

■  :  equation  117. 

value  cdileu  by  LKAf  to  factor  die  material 
into  an  upper  triangular  matrix  described  in 


-  .  a  .ibrnut  ins  called  by  C PLATE  to  determine  the  local  to 

'  ..  :a,  *  1  :o--mation  matrix  as  described  by  equations  76 

.  The  routine  checks  for  proper  area  definitions 

-  :  s  i;t ire  called  uy  hOCAL  to  perform  vector  cross 

:  ir,  equations  76,  78  and  79.  The  resulting 

t  '  r  compont/r: :  a  are  normalised  to  unit  vectors. 

-  ..  /coroutine  Called  by  LOCAL  to  perform  an  in-place 
:  c..  re  matrix  transpose  of  the  local  to  global  transfor¬ 
mation  '  oe^caon  c 0. 

u:c  trirx  .lied  by  LOCAL,  to  transform  global  element 

.  ■  ;  :  i  •:  *  :  :  >  in  c  a  1  ele  me  n  l  o  o  o  r  d  inat^s.  Transform  a  i  1  o  n 

-i :  i  i  x  .  ---  w.- .  i  by  e  q  eat  ion  80. 

..  s  aiS'O.i r;  i.n<  anile:  Lv  CLLATL  to  ^et  matrix  array  space  to 
•tnv  vj.ua.  L;.u-o ;  :  i  ca1*  1  y ,  it  in  a  red  to  set  matrix  space 

i  .  ■  . '  e  r  *  o  • 

s  -  -.4  . r  >ut  i  no  called  by  C PLATT  to  form  composite  plait; 

.  ic,.,r  ;/  train  a  i.  tola  'emerv?  :  unci  1  on  A  as  described  bv 

*  «.  5  i 

•  nation  n?  .i  ;;'l  anon  equation:,  *'  tr.^'Uyh  A  5 . 

-  .  a  :  *!  in-  called  bv  CLuATT  to  ,u  e  jnalatr  1  :»e  rv:,Ui  is 
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-  -•  -  ■  calico  oy  LjaLRN  to  porrorm  a  special 

.rr  r  ix  ::*u  - r  .  pi  y  where  the  leading  matrix  is  an  upper 

aiTix ;  descritcj  oy  equation  v b . 

;  i.ubrout  i  no  .  a’  ioc.  by  UbKERE  to  perform  a  special  sym- 
■tdic  multiply  .1:--  aeeooa  in  equation  96. 

-  1  -a; '.‘O,.:  in  a  eulleo  by  /PLATE  to  perform  a  lower  sym- 

Cine  te.P.  . multiply  o:  the  natural  stiffness  components 
-nrirv  'i«  ■  n.orio.xi  integration,  as  in  equation  108. 

-  a  0  u!  r\-.d  .  s  ■  cm.:  led  by  Cl  LATE  to  perrorm  the  lower 

vim:.*  ^  .  •  r  ..melon  o:  tne  stiffness  matrix  components 
uritr1  i...  ...  .0.1  integration ,  as  described  by  equation  ICS, 

-  a  a  a:  :  o .. :  cal  lea  by  C  PLATE  to  determine  the  distri¬ 

bution  c;  u:ea,  bp  node,  for  the  quadrilateral  plate  elc- 
=en..  In.la  area  function  is  needed  in  equations  112  and  11 

-  a  a  .fr-nou ;  ine  called  by  CP  LATE  to  form  a  transformation 
.atnix  f r  mr  nrPa ra^  to  local  coordinates  as  described  in 


-  ;usr  *'>  :t  ine  called  by  CP  LATE  to  form  the  natural  to 
iota!  coord  1 nate  1  ran s format ion  as  shown  in  equation  67. 
t.  o  t  rune  formation  matrix  is  diagonal,  only  the 
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g  condit ions . 

was  compared  with  many  simple  degenerate 
e ems  and  was  found  to  produce  excelled 
were  performed  for  element  TYPE  1  an 
were  more  favorable  for  the  simple  cases, 
•.suit.;  in  larger  displacements  while  TYPE 
.ewer  than  exact  solution  values, 
ret  tier,  describes  a  group  of  verification 
.nrvec  shells  and  a  doubly  curved  blade. 
tv  obtained  by  classical  plate  and  snell 
:  tiUtc  mater!  .1 1  properties..  The  pirn-,  eburt 
1. 1 ■  r  see  i.e s  t  o  a p p r o x  t  m a t e  . . < > .i  a t  x  o : . s  v_ - r  v  a r 
v  so  mi t sons .  Once  tne  A,  o  mat r  i c.e s  w 
soy  c  o:\k1i  V :  ens  applied,  an  approx  -  mate  s: 

.1  never  or  s  is. ht—  term  ^xr  ts.^cn .  Most 
o;s.r  slonentn  per  is  nstm.  ronv^ryo's 
■  ■  :,ts i od  :.s  nc e  a  g ; . «- ; - <  *  . . :  u o y  > ;■  n  t r. t  e ^ <... 

;  erenoe  I  . 


r ‘.n/Oo/so]  under  uniformly  distributed  lateral  load 
with  all  edges  simply-supported 


.size  of  la;ninated  plate: 


Deflection  at  center,  W 


(  A  ) 


( W  ) 


62.38  x  10 
62.34  x  10 


c; 

-4 

-4 


Po 


3 

/ib 

3 

/lb 


exact 


4.  L  Go/ 902  Jt  UNDER  UNIFORMLY  DISTRIBUTED  LATERAL  LOAD 
WITH  ALL  EDGES  SIMPLY-SUPPORTED 

(i)  Size  of  laminated  plate: 

a  =  10"  b  =  10"  t  =  4h  =  .2" 


(2)  Properties  of  plate: 


.  0275 

lb..  3 
/in 

’3.  3 

#  2 

,  2 

3.3 

0 

0 

“.135 

0 

0 

-.135 

0 

0 

r 11 

.67 

.  67 

11 

)6  lb/in 


1  Ox  10J  lb-in 
0  .  67  I 


i,r, 


(3)  Loading  condition: 
p(x,y)  =  p^ 


FIGURE  11 


(4)  Boundary  conditions: 

x  =  0 ,  a ;  w=0  M  =0  v 

x 


0  N 


0,  b;  w  =  0  M  =0  u  =  0  N 


(5)  Deflection  at  center,  Wc ; 

(W  )  =  114.4  x  10“4  p  in3/lb. 

c  exact  *o 

(W  )  =  113.6  X  10-4  d  in  /lb. 


o.  [G/30]t  UNDER  IN-PLANE  LOAD  WITH  TWO  EDGES  PERPENDICULAR 
TO  THE  DIRECTION  OF  LOAD  FREE  AND  OTHER  TWO  EDGES  SIMPLY 

SUPPORTED 


C) 


o  =  1  psi 

X 

(4)  Boundary  conditions: 
x  =  0  ;  u  =  0 
x  =  a  ;  u  i  0 


(5)  Curvature  at  x-direction, 

(K  )  =  . 3 0 5o 

x  exact  x 

(K  )  =  -.30570 

x  sap.  x 


[SO/0-?/90]-r  CURVED  PLATE  UNDER  UNIFORM  PRESSURE  WITH  ALL 
EDGES'  SIMP L.Y -  SUP PO RTED 


M  cn 


j.  NATURAL  FREQUENCY  OF  A  TYPICAL  BLADE  CONFIGURATION 

(1)  Type  of  Blade:  (Figure  16) 

A  typical  blade  used  in  aircraft  engines,  named  J79  B/AL, 
was  approximated  using  64  quadrilateral  plate  elements  in  SAP4A 
(both  TYPE  6  and  TYPE  9).  An  equivalent  NASTRAN  model  was  also 
run  along  with  an  experimental  test  to  determine  fundamental 
frequency  at  zero  frequence  (results  from  AFSC  -  Wright-Patterson 
AFB)  . 

(2)  Properties  of  blade: 

Material  used  corresponded  to  the  input  used  in  the  NASTRAN 
run  using  anisotropic  material  properties. 

Leading  edge: 

r  =  .000407  lb  sec^/in  9 

Cxx  =  26.9E8  ib/in2  Cxy  =  4.6E7  lb/in^ 

Oyy  =  20.3E8  lb/in2  Gxy  =  7.3E7  lb/in2 

Blade: 

2  . 

p  =  .000261  lb  sec  /in 

same  material  coefficients  as  above. 

(3)  Loading  Condition: 

Mass  and  Stiffness  distributions  for  eigensolut ion . 

(4)  Boundary  Conditions: 

Base  of  blade  completely  fixed  and  rest  of  blade  free. 

(9)  Natural  Frequency  of  Blade:  (Fiest  FLEX) 

(f)exp  =  110  Hz. 

( 1  )Nastr..n  r  i  06  ,  7  Hz  . 

(f )sap  (TO )  =  i08 . 3  Hz. 

(f)sap  (T9)  =  112.8  Hz. 
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SECTION  V 

DISCUSSION  AND  CONCLUSIONS 


A  quadri lateral  composite  plate  finite  element  has  been 
a -idea  to  the  SAP  IV  computer  program  library  to  be  used  on  plate 
tvpe  structures.  The  element  is  a  four-noded  sub-parametric 
;  lat*  plate  element  excluding  transverse  shear  deformations.  The 
clement  is  "incompatible"  relative  to  mid-plane  surface  rota¬ 
tions  along  the  inter-element  boundaries.  The  element  converges 
relatively  well  despite  the  incompatibility,  as  long  as  the 
element  maintains  a  relatively  rectangular  shape  with  a  reasonable 
element  aspect  ratio. 

The  element  was  run,  modelling  various  simple  plate/beam 
configurations  and  showed  excellent  results.  More  complex  models 
were  designed  to  test  the  composite  laminate  behavior  of  the 
clement,  as  described  in  section  IV.  The  models  included  flat, 
cylindrical  and  doubly-curved  shells.  The  results  obtained 
compared  favorably  with  classical  series  solutions,  (approximately 
1-6%  disagreement)  The  results  obtained  in  section  4.9  for  the 
typical  blade  conf igurat ion  shows  that  the  NASTRAN  and  SAP 
(TYPE  6)  element  values  were  slightly  lower  than  the  experimental 
number  while  the  result  of  the  new  element  (TYPE  9)  was  slightly 
higher.  The  new  element  is  based  on  an  incompatible  formulation 
and  in  general,  convergence  is  not  guaranteed.  But,  in  more 
oases,  the  element  has  been  found  to  be  slightly  stiffer  than 
compatible  elements  and  therefore  the  frequency  is  higher. 

The  quadrilateral  element  will  be  inherently  stiff  if  the 
four  points  do  not  represent  a  flat  surface.  Therefore  the  ’‘quad” 
element  was  relaxed  to  better  represent  shell  behavior  by  allowing, 
the  flat  plate  element  stiffness  coefficients  to  be  transformed 
relative  to  the  local  nodal  coordinates  of  the  element. 

The  thin  plate  theory  used  to  develop  the  element  does  not 
account  lor  normal  torsional  effects.  Therelore,  non  global 
adjacent  elements,  when  assembled,  will  produce  a  singularity 


normal  to  the  plate  if  the  two  elements  are  coplanar.  To  avoid 
this  singularity  in  general,  an  artificial  torsional  stiffness 
or  scaffolding  matrix  was  added  to  each  of  the  four  nodes. 

This  does  violate  element  equilibrium  but,  if  the  magnitude  of 
the  coefficients  are  maintained  relatively  "soft"  compared  to 
the  plate  bending  characteristics,  overall  equilibrium  is  closely 
m.i  -  nta  ined  . 

The  composite  plate  element  can  be  used  in  the  existing 
static  and  dynamic  analyses  contained  within  the  SAP  TV  program . 
it  cun  be  used  with  all  the  existing  elements  in  the  finite  ele¬ 
ment  library  as  long  as  the  model  effects  are  correct. 

The  element  was  not  being  developed  to  degenerate  to  a 
triangular  plate  element.  If  the  element  is  used  as  triangular, 
the  local  effect  is  ’’too  stiff" .  If  the  fourth  of  the  quadri¬ 
lateral  nodes  is  placed  mid-plane  on  a  triangular  side,  a  better 
approximation  can  be  obtained. 

The  composite  element  presently  does  not  contain  geometric 
tillering  effects  ;>uch  as  those  required  in  high  speed  centri- 
:  u-v  machinery. 

further  Developments : 

further  work  ;  hould  be  done  to  include  the  geometric  st i : - 
inning  c  'Vo  :  :  ie  1  cm  ■  to  account  for  centripetal  acceleration 
♦  ■ ;  feet  s  o :  h i  g h  s pin  i >  1  a d e  s y s t em s  .  Th i  s  go ome t r i c  matrix  of 
coefficients  would  allow  a  better  approximation  of  blade  bending 
stresses  and  closet*  representat ion  of  blade  natural  frequencio: 
at  high  spin. 

A  p re- processor  program  should  be  developed  to  handle  complex 
material  laminates  as  a  function  of  blade  position.  This  pro¬ 
gram  would  compute  the  A,  B  and  D  matrices  and  thermal  vector: 
needed  in  the  SAP  program  and  produce  a  storage  file  for  stress 
recovery.  This  program  should  also  plot  all  information  i  or 
input  chec k  Lng . 


A  pow  should  Le  written  to  retrieve  stress  in- 

ion  and  material  information,  by  lamina,  to  be  used  with 
arm. it  ion  output  and  curvatures  at  the  mid-plane  surfaces  to 
■;\duce  individual  lamina  stress  to  be  used  in  failure  criteria. 
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appendices 


AsrEE.UX  A  -  1  ii PUT  TO  ELEMENT  TYPE  9  IP  SAr4  A 

A.«.  comp] o  t e  input  to  the  SAP  IV  program  ;  n  descrineu  i  n 
v.: .  u’tna1  I.  Tnc  :n-w  element  AYAj;  9  defined  in  PAi  9A  a 


: udcd  ana  snoulJ  De  appended  to  reference  1.  The  in  m 

•'  aiv.rnt::’.’:  TYPE  9  in  similar  to  element  TYPE  6.  In  :  act ,  o.« 


aixn  replace  TYPE 


. ne  rormat  o:  t..e  .in  .1 


t  net  .:,  tion  is  consistent  with  the  SAP  IV  input.  Varieties 
.Yelled  integer  must  be  right  justified  and  floating  po-r.t 
ari abies  should  include  a  decimal. 


"Mi  -  ITE  PLATE  INPUT  TO  SAP4A 


TYPE  9  Composite  Plate  Element  (QUADRILATERAL) 

Note  Columns  Variable  Remark 

A.  Control  Card  (615,  110) 

5  NPAR(l)  Number  S 

6-1C  NPAR(2)  Number  of  Plate  Elements 

11-15  NPAR(3)  Number  of  different  materials 

(1)  16-20  NPAR ( 4 )  Material  Type  Key 

=0  Composite  Material  Prop. 

=1  Standard  Anisotropic  prop¬ 
erties  (same  as  TYPE  6) 

(2)  21-25  NPAR(5)  Number  of  Global  Material  Vectors 

If  zero  or  blank  then  global  X 
direction  is  assumed  to  be 
material  x  axis. 

26-30  NPAR(6)  Integration  Order  (default  set 

to  2 ) 

(3)  31-40  NPAR(7)  Rotation  Stiffness  Factor 

(integer  number) 

B.  Material  Property  Information 

Two  types  of  material  can  be 
input  to  element  type  9:  general 
composite  material  and  aniso¬ 
tropic  material. 

o.i  Composite  Material  Properties  (NPAR (4 ) .EQ. 0) 

Five  cards  must  be  input  for 
every  different  material.  (NPAR(3)) 


Card  1 : 

(110, 

20X  , 

4F1Q. 

0) 

:  -io 

NN 

Material  identification  number 

11-30 

Blank 

31-40 

DEN 

Mass  density 

41-50 

AT  ( 1 ) 

51-60 

61-70 

AT  ( 2  ) 
AT(  3 ) 

A„ 

thermal  vector 
(equation  66) 

components 

Card  2  : 

( 6F10 . 

0) 

1-10 

BT(  1 ) 

11-20 

BT  ( 2  ) 

21-30 

BT  (  3  ) 

B„ 

-  i 

thermal  vector 

component : . 

31-40 

DT(  1 ) 

(equation  67) 

41-50 

DT  ( 2  ) 

thermal  vector 

components 

51-60 

DT(  3  ) 

—  1 

(equation  68) 
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3  : 


(6F10.0) 


:  - 1  o 

A  (  1  ,  1 ) 

11-20 

A( 1 , 2 ) 

A 

21-30 

A(l,  3) 

3 1  -  4  C 

AC 2 ,2) 

41-50 

AC  2  ,  3  ) 

51-60 

A(3  ,3) 

Card  4 ; 

(6F10.0) 

1-1C 

B  ( 1 , 1 ) 

11-20 

B(  1 , 2  ) 

B 

21-30 

B  ( 1 , 3  ) 

**' 

31-40 

B(  2 , 2  ) 

41-50 

B(  2 , 3  ) 

51-60 

B  (  3 , 3  ) 

Card  5 : 

(6F10.0) 

1-10 

D(l,l) 

11-20 

D  ( 1 , 2  ) 

D 

21-30 

D(1 , 3 ) 

31-40 

DC2 ,2) 

41-50 

DC  2 , 3  ) 

51-60 

DC  3 , 3  ) 

Matrix  coefficients  (upper 
triangular)  (equation  63) 


Matrix  coefficients  (upper 
triangular)  (equation  64) 


Matrix  coefficients  (upper 
triangular)  (equation  65) 


*2  Anisotropic  Material  Properties  (NPARC 4 ) . EQ . 1 ) 

4)  Two  cards  must  be  input  for  every 

different  material  (NPAR(3)) 


Card  1 : 

(110,  2  OX , 

4F10.0) 

1-10 

NN 

Material  identification  number 

11-30 

Blank 

31-40 

DEN 

Mass  density 

41-50 

AX 

Thermal  expansion  coefficient  ax 

51-60 

AY 

Thermal  expansion  coefficient  ay 

61-70 

AXY 

Thermal  expansion  coefficient  axy 

Card  2  : 

(6F10.0) 

1-10 

CXX 

Elasticity  element  Cxx 

11-23 

CXY 

Elasticity  element  Cxy 

21-30 

CXS 

Elasticity  element  Cxs 

31-40 

CYY 

Elasticity  element  Cyy 

41-50 

CYS 

Elasticity  element  Cys 

51-60 

GXY 

Elasticity  element  Gxy 
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C.  Global  Material  Vectors  (25,  5X ,  3F10.0) 


NPAR(5)  material  vectors  must  be  input  (except  if 
zero  or  blank) 


1-5 

NV 

Material  vector  identification 
number 

6-10 

Blank 

11-20 

DX 

X  direction  cosine 

21-30 

DY 

Y  direction  cosine 

31-^0 

DZ 

Z  direction  cosine 

it  Load 

Multipliers  (5 

cards ) 

Card  1 : 

(4F10.0) 

1-10 

PA 

Distributed  lateral  load  multi¬ 
plier  for  load  case  A 

11-20 

PS 

Distributed  lateral  load  multi¬ 
plier  for  load  case  B 

21-30 

PC 

Distributed  lateral  load  multi¬ 
plier  for  load  case  C 

3 1  -  4  C 

r  D 

Distributed  lateral  load  multi¬ 
plier  for  load  case  D 

Card  2 : 

:  (4F10.0) 

1-10 

TA 

Temperature  multiplier  for  load 
case  A 

11-20 

TB 

Temperature  multiplier  for  load 
case  B 

21-30 

TC 

Temperature  multiplier  for  load 
case  C 

31-40 

TD 

Temperature  multiplier  for  load 
case  D 

Card  3: 

:  (4F10.0) 

1-10 

XA 

X-direction  acceleration  for 
load  case  A 

11-20 

XB 

X-direction  acceleration  for 
load  case  B 

21-30 

xc 

X-direction  acceleration  for 
load  case  C 

31-40 

XD 

X-direction  acceleration  for 
load  case  D 
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Card  : 

( *  i  F  i'j.  0  ) 

1-10 

A 

Y-di rect ion 
load  case  A 

acce  1  t-rat  ion 

1 

fO 

o 

YB 

Y-c ‘ rect ion 

acce 1 erat ion 

Loud  C.iuO  p 

21-30 

YC 

Y-d i rect i on 

acce  1  t  ion 

lead  care  0 

31-40 

V  • 

Y-d  i  rect  1  or: 

acce lerat i on 

load  case  D 


Card  5  : 

(  4  f  i  G  .  0  ) 

1-10 

ZA 

Z-direct ion 
load  case  A 

acceleration  for 

11-20 

ZB 

Z-direct ion 
load  case  B 

acceleration  for 

21-30 

ZC 

Z-direction 
load  case  C 

acceleration  :  or 

31-40 

ZD 

Z-direction 
load  case  D 

acceleration  f or 

I! lament  Cards  CITS,  12,  13,  12,  13,  15,  4F10.0) 
One  card  for  each  Ni:‘Ak(2)  clement . 


1-5 

NN 

Element 

number 

(5) 

0  -  i  0 

I 

Node  I 

11-1  b 

J 

Node  iT 

lb- 2  0 

K 

Node  K 

21-25 

L 

Node  E 

4  l;  ) 

2  G  -  2  7 

NO 

No.  of 

u a  i  i  s s  integration  :  c !  n  t  . 

(7) 

2  8  -  3  0 

IV 

Mater ial  vector  ident  i  f  i  cat  • 

number 

(d) 

3  1  -  3  2 

ir-.iisi: 

Previou 

s  Clement  re-use  code 

"  0  new 

e lement 

(9) 

33-35 

IM 

-1  use  previous  element 
Material  identification  number 

(13) 
(11  ) 

36-40 

41-50 

INCL 

TH 

Clement 

Clement 

generation  parameter 
thickness 

(12) 

51-60 

PR 

Clement 

lateral  pressure 

(13) 

61-70 

TO 

Mean  temperature  variation  from 

the  ref 

erence  level  in  uncie- 

t ormed 

pos i t i on . 

71-80 

TO 

Mean  to 

m]  >e  ra  t  u r e  g r a d  i  c  n 1  acros 

the  she 

11  tfii  ekne  ss  . 

Ndtr; 

. : 

U  ) 

ill  cine  nt  TYPE 

4  allows  two 

•  lit  I  e re n  1 

i  t  >r*Tii  ■  d  input.  The 

:  irst  form  is  tor  laminate  matrices  whi!<  the  .Aicomi  i*.  t  oj 
anisotropic  matrices.  The  Later  form  is  identical  to  ele¬ 
ment;  TYPE  6  input. 


hi 


(2)  A  material  plobal  axis  must  be  defined  relative  to  the 
material  properties  formed  in  the  A,  B  and  D  matrices. 

(3)  Kt.»t.at  iona  1  Stillness  Factor  is  set  by  multiplying  NFAR(V) 
times  l.F-8.  Default  for  N  i  Al<  ( 7  )  is  100. 

(4)  Material  input  in  this  section  is  identical  to  that  of 
element  TYPL  f* . 

(5)  The  T,J,K  and  I,  indices  detine  the  element  connectivity 
and  also  the  element  normal.  The  element  "z"  coordinate 
is  formed  bv  the  ripht  hand  rule  as  I  poes  to  J  p  oes  to  K , 
etc.  The  element  local  axis*  is  determined  by  the  projec¬ 
tion  of  the  plobaL  material  axis  onto  the  element.  Once 
the  x-ax is  is  determined,  the  local  y  is  formed  from  z 

and  x.  All  stress  output  is  in  this  relerence.  Tf  node  L 
equal  K  or  i!  zero  or  left  blank,  the  program  will  assume 
that  the  element  is  tri  anpu  1  ar .  The  result  inp.  local  stiff¬ 
ness  is  then  "too  stiff". 

(<  )  The  number  of  Mauss  integration  points  can  vary  as  2  or  3. 

It  a  value  is  set  above  or  below  these  numbers,  the  prop, ram 
will  reset  it  to  2.  The  default  value?  i :  ^  set  NPAR(6). 

(7)  The  ploLvil  material  vet' tor1  must  be  p, renter  than  or  equal 

to  1  and  1  ess  than  or  equal  to  NPAK(S).  It  NPAR(S)  is  blank 
or  zero,  then  N  PARC  5)  is  set  equal  to  1  and  the  plobal  vec¬ 
tor'  is.  alipned  alonp.  the  plobal  X  axis.  Default  value  if; 
set  to  1 . 

(8)  if  an  element  has.  the  same  planar  size,  lame  oriental  ion  in 
space  and  the  same  element  load  inp,  parameter:;  as  the  pre¬ 
vious  e.le;:ient  ,  then  sett  inp  IKFWSI!  equal  to  1  wil  l  use  the 
same  pi  o  b  1 1  element  stiffness*  and  load  vectors;  for  assembly. 
Default  is  set  to  0. 

(3)  The  mater' i a. 1  ID  number  must  be  between  1  and  NPAR(3).  De¬ 
fault  is  set  to  1. 

(ID)  clement  Meneration  Parameter:  blcrmnt  cards  must  be  in 

'■lement  number  sequeno<  * .  It  element  car*ds.  are  omitted,  the 
prop  ram  will  penepate  the*  miss  inp,  rcirds  as  follow:;: 

The  increment  for  the  element  number  is  one. 

i.e.,  NM.  _  -  NfJ .  +  1 

Th'*  f  ■  f  i  r  ‘fe  ■ ,  p< » r  e  I  i  n  s  i  iirniii'iit  l*»r  nodal  r^nnert  ivity  i:.  INTI.. 


■’  p  *** 

-i+1 

= 

I. 

i 

INCL 

T 

"  i  +1 

= 

J. 

i 

+ 

INCL 

Ki+i 

= 

K. 

l 

+ 

INCL 

Li+i 

= 

L. 

l 

+ 

INCL 

If  INCL  is  left  blank 

then 

INCL  is 

set  to  1 

. 

Material 

identification,  element 

thi 

ckness , 

distributed 

latera I 

load ,  temperature 

and 

temperature 

gradient 

for 

the  element 

are  then  the  same 

as 

for 

th 

e  first 

element 

in 

the  generated 

group.  The  last  element  card  must  be  input  to  exit  element 
group  properly. 

(11)  The  plate  thickness  is  used  mainly  to  compute  the  mass  of 
the  element.  Default  is  set  to  1.0. 

(12)  The  pressure  is  normal  to  the  surface  of  the  element.  The 
positive  direction  for  the  pressure  loading  vector  is  in 
the  positive  direction  of  the  local  z  coordinate. 

(13)  The  temperature  required  is  the  mean  temperature  difference 
(TO)  from  the  reference  temperature  of  the  element  in  a 
undeformed  state.  The  TG  is  the  mean  thermal  gradient 
through  the  element  thickness. 


APPENDIX  B 


INPUT  TO  PRL-PROCL’SSO R  PROGRAM  LAYUP 


A  iT'e-proct ■  :v.;or  program  waf;  developed  to  calculate  the 
,  B  and  D  matrices  in  described  in  equations  6  3  through  6  3  and 
e  r  hernial  load  vi\.  i  ora  described  in  equ(«t  ion:,  66  throu^n  68. 
:ie  material  Lnt ormat ion  by  layer  in  input  no  lamina  position 
nd  : iber  nriuntat ion .  Tne  pre-processor  then  computes  the  0 
it  rices  invoivin.  the  tensor  transformation  and  constructs  the 
loment  materia1,  matrix  C  relative  to  the  mid-plane  of  the 
late.  Similar  thermal  vectors  are  also  computed  and  output, 
his  information  is  then  used  directly  for  SAP4A  -  TYPE  9. 


INPUT  TO  PROGRAM  LAYUP 


A.  Number  of  Cases  Card  (12) 

Note  Columns  Variable  Remark 


(1) 

1-2 

NCASES 

5 . 

Heading  or  Tit 

le  Card  (8A10) 

(2) 

1-60 

ITITLE 

C. 

Laminate  Data 

Card  (415) 

(3  ) 

l-o 

NLAM 

(4) 

o  - 10 

NMAT 

11-15 

ITRAN 

16-20 

IORD 

.  > 

Material  Prope 

rty  Card  (6F10 

(  b  ) 

1-10 

EI1(K) 

■’  i  _  ■■)  ,'j 

j:22(K) 

;  :  -  j  0 

XN1(K) 

3.-40 

012 (K) 

S  I  -  5  0 

THERMl(K) 

51-60 

THERM 2 ( K ) 

Enter  the  total  number  of 
laminate  configurations  to  be 
considered . 


Enter  the  title  information  to 
be  printed  with  the  output. 


Enter  the  number  of  laminae  in 
the  laminate. 

Enter  the  number  of  different 
materials  in  the  laminate. 

If  ITRAN  is  zero  or  blank  the 
transformed  thermal  properties 
of  each  lamina  will  not  be  part 
of  the  output. 

If  IORD  is  left  blank  the  lami¬ 
nate  ordinates  will  not  be  output 

) 

Enter  lamina  modulus  in  fiber 
direction 

Enter  lamina  modulus,  in  direction 
transverse  to  fibers. 

Enter  manor  Poisson's  ratio. 

Enter  lamina  shear  modulus. 

Enter  C.T.l.  In  liber*  direction. 
Enter  C.T.E.  in  transverse 
direction . 


L . 


Lamina  Data  Card  (P10.3,  110,  FI 0.3) 


1-1C  T(J) 

11-20  MATH  J ) 


Enter  thickness  of  the  J  lamina 
Enter  the  number  that t^dent if ies 
the  material  of  the  ‘  lamina.^ 
Enter  the  orientation  of  the  J  il 
lamina  with  respect  to  the 
lamina  t  e  axe:.; . 


21-30 


PJI(J) 


Notes : 


(1)  There  are  no  program  restrictions  on  the  number  of  cases 
that  may  be  analyzed  in  a  single  run. 

(2)  Begin  each  new  data  case  with  a  heading  card. 

(3)  The  program  is  currently  capable  of  handling  up  to  48  lami¬ 
nae  per  layup.  This  can  be  increased  by  changing  the 
appropriate  dimension  statements  as  shown  in  the  program 
listing. 

(4)  Although  no  limitation  on  the  number  of  different  materials 
need  be  imposed,  the  program  dimension  statements  currently 
allow  for  a  maximum  of  NMAT  =  4.  However,  this  can  also  be 
increased  if  necessary. 

(5)  This  material  data  is  vendor  information.  One  card  is 
required  for  each  material  (i.e.,  K  =  1,  NMAT). 

(6)  One  card  is  required,  per  layer,  in  the  laminate.  (i.e., 

J  =  1,  NLAM) . 
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APPENDIX  C 


COMPUTER  RUN  -  INPUT  AND  OUTPUT 


This  appendix  contains  the  complete  input  and  resulting 
output  of  the  SAP4A  program  using  the  new  composite  plate  ele¬ 
ment  (TYPE  9).  The  example  is  the  model  contained  in  section 
4.7  (Curved  Plate  under  Uniform  Pressure).  The  first  section 
contains  a  listing  of  the  card  images  used  to  execute  the  pro¬ 
gram.  The  second  section  is  the  complete  output  resulting 
from  this  input. 
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O  O  O  O  O  O  O  O  O  O  G  O  O  O  CD  O  O  ' 
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00000000000000(3  00  0 


X  o  lO  ro  x  x  O  un  oo  x  x  O  U"i  no  x  x  O  zv  O 
go  OLO^Lnaio^«d-ir)^oif)^^aOLn^ 
uj  ocvj^LOir>oc\j<tLOLno(\i^inir50CMC 


d  o  o  o  o 


Q 

cl 

o 

0>-00000000000000  0000 
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CD 

o 

o' 

GO 

cn 

UJ 

\ 

z 

CD 

O 

LO 

r — 

o  o 

o 

0 

0 

O 

o 

X 

I — 

0 

0 

o 

1 — 

0 

0 

0 

O 

1 — 

0 

0 

o 

o 

i - 

o 

0 

0 

cr\ 

o 

cm 

CD 

X 

1 

v 

Z 

o 

►—4 

S 

O 

Z 

- _ - 

it 

II 

it 

II 

II 

z 

ZD 

► — < 

II 

II 

ll 

II 

o 

t— 

O  LU 

CL 

/ - -X 

L-t 

X 

1 — 

o 

o 

o 

o 

1 

0 

0 

o 

o 

i — 

o 

o 

o 

o 

, — 

o 

o 

; 

*— »  GO 

1 — 

d 

- - N. 

:> 

!— 

X 

p. — . 

d 

GO 

h“  Z 

O 

CL 

o 

GO 

►— t 

d 

GO 

LU 

o  o 

LU 

O 

d 

o 

d 

Q 

_ 1 

t— 

CL 

GO 

d  CL 

CL 

LU 

. — >x 

z 

I - 

1 — 

Z 

cl 

Z 

X 

GO 

LU 

- — - 

Od  GO 

GO 

t — 

X 

z 

d 

o 

X 

F— 

0 

0 

0 

0 

0 

0 

l 

o 

o 

i — 

0 

0 

o 

o 

r - 

o 

o 

cl 

*— < 

I — 

LU 

t — i 

z 

h—  LU 

Z 

LU 

z  o 

— _ ^ 

O 

CD 

1 

o 

o 

GO 

O 

x 

x  or 

LU 

► — i 

O  o  LU 

LU  GO 

CD 

LU 

o 

CL 

h-  < 

z 

o 

LU 

GO 

o 

O  Od 

o 

LU 

1— 

X 

> 

z 

o 

LU 

z 

o  o 

Z 

h- 

51 

1—  O  d  o 

3 

CD 

3 

CL 

X 

1 — 

o 

0 

0 

0 

0 

0 

l 

0 

0 

r— 

0 

0 

o 

o 

r— 

o 

o 

a: 

u_ 

_ 1 

LU 

o 

_ - 

M  U  LU 

O 

O 

■ — «• 

ZD 

CL  I— 

CQ 

d 

Q_ 

d 

o 

d 

2:00^ 

h-  d  o 

CL 

LU 

o  d 

GO  O 

3 

Z 

o 

o 

Z 

Q 

LU  d  LU 

LU 

dao: 

cn 

oc 

LU 

LU  1— 

CQ  LU 

CL 

LU 

z 

o 

_ 1 

O  CL 

a 

1—0  0 

LU 

k— < 

o 

x  d 

3  > 

LU 

3 

X 

F - 

0 

0 

o 

1 — 

0 

0 

0 

o 

r— 

o 

o 

o 

o 

, — 

o 

o 

o 

LU 

»— t 

z 

LU 

_ 1 

Ll 

O  GO  2l  Ll 

Od 

O  O 

LU  O 

GO 

CL 

LU 

1 — 

o 

1 

GO 

O 

s: 

Z 

> 

z 

CQ 

< 

Ll 

Lx- 

Ll 

Ll 

U-  O 

GO 

d 

>— 1 

O 

_l 

o  o  o  o 

GO 

^  ^  • 

rv 

z 

* 

o 

z 

GO 

o 

t— t 

O' —  CVl 

ro 

d 

o 

o  — 

1 — 

o 

CL 

i-  — 

o 

cl 

or 

or 

CL 

GO 

•  •  . 

. 

.  • 

CL  < 

►— t 

o 

Od 

gD 

LU 

LU 

LU 

LU 

>-  aaa 

cr 

O  1— 

O'  O'  LU  OC 

h- 

F— 

_ 1 

LU 

LU 

a: 

CQ 

ao 

CQ 

CQ 

_j 

LU  LU  LU 

LU 

LU 

3 

LU  LU 

CQ  LU 

d 

LU 

d 

LU 

CQ 

1 - 

CO 

d 

LT) 

6 

7 

8 

CTi 

o 

f— 

X 

CO 

d 

LiO 

lO 

X 

cc 

H- 

s: 

s: 

s: 

d 

s:  i— 

3 

CL 

o 

o 

5: 

r— 

f— 

r— 

r— 

r— 

f— 

r— 

o 

o 

3) 

3 

o 

z 

o 

3  Z 

O' 

d 

o 

o 

3 

d 

z 

z 

Z 

z 

d 

GO 

z  •— « 

LU 

1— 

z 

z 

Z 

d 

z 

CL 

o 

d 

GO 

o 

71 


O  O  O  O  O  O  O  OOOOOOOOOOOOOOOOOOOOOOOOO 

ooooooo  ooooooooooooooooooooooooo 

ooooooo  ooooooooooooooooooooooooo 


ooooooo  ooooooooooooooooooooooooo 


ooooooo 


rx  c-*  o  o  oo  r>-  n 

if)  CTi  O  Lf)  ^  lD  (Ti 
LD  ld  o  cv  ^  ID  Ul 


NOLnrONNO^CONNOLOnNNOlDnNNCLOOONN 

uo  OLn^^^o^^L^c^oLn^Lnaoin^in^o^c^cr. 

lu  ocvj^Ln^oovj^LnLDocM^LDLDOCM^Lmnoosj^Lnir; 


o 

a: 

o 

OOOOOOO  0>-0000000000000000000000000 

OuDOOOOO  O  OOOOOtD-OLOLOLnoOOOO^OUOLnLOLOOOOOO 

X  C-  O  O  O  O  O  OOOOOCMCMWWCMLT)Lf)Ln^^NNNNNOOOCO 


n  lTuo  id  lo  o  ooooo* —  i —  » —  « —  r-c\jc\jcMC\jc\irororoooroL/)ir)ir)i^LO 

O 

o 

OOOOOOO  —JXOOOOOOOOOOOOOOOOOOOOOOOOO 

looouoo^-oo  c  OLDOLOOOLnouoooLOoi-noOLOOunoOLOOuOO 

ix  o  o  cm  lo  r--.  o  o  ocMLnr^oocMLOixoocMLnrxoocNjLOrxoocMLOxo 

.  o  . 

CO  LD  o  » —  CM  CO  1-0  O  O  < —  C\J  CO  o  o  < —  CM  CO  LO  O  < —  CM  CO  UO  O  i —  CMOOLDO' —  CM  CO  O 


o  <—  o  o  o  o*— 

O'—OOOOr— 

oo-ooo^- 

J 

O  o  r-  o  o  o  o 

O  O  «—  O  O  O  f— 

I 

o  —  o  o  o  o  •— 

cr>  o  » —  cm  co  lo 

I —  CM  CM  CM  CM  CM  CM 


r-v,  O  O  O  O  < —  OOOOr—  OOOOl —  OOOOl —  o  O  O  O  I — 
CM 

uo 

LlJ 

o 

O  1 —  1  '  *  1 —  OOOOr—  OOOOl —  OOOOr—  oooo* — 

O  >-  l  I  \  1 


OOOOf—OOOO' —  oooo< —  oooo< —  < —  • —  t —  i — 

I  I  I  I 


«—  •—  O  O  O  O  *—  OOOOr—  oooo^oooo 
I  I  I  I 


0000^-0000*—  OOOO*—  OOOO*—  »—  r—  r- 

till 


OOOO*—  oooo*— oooo*— 

I  I  I 


H— 

C  oc 

cn  UJ 

UJ  UJ  CO  *—  CMCO^TLOOJXCOOO 

z  o  s:  <— 

UJ  o  o 

O  O  2T 


Oon^iniDNCOWO' —  CMCO^fLO 
■ —  I —  I —  I —  1 —  I —  f —  I —  CM  CM  CM  CM  CM  CM 


72 


EQUATION  NUMBERS 


x 

_ _ 

x 

o 

o 

LO 

. —  x 

co 

o 

cn 

if)  t —  rx 

o 

CO 

cn 

LO 

, — 

o 

X 

r— 

LO 

cn 

o 

x 

CXI  CM 

CO 

CO 

cn  in 

LO 

LO 

X  CO 

00 

cn 

cn 

cn 

>- 

o 

o 

o 

o 

o 

o  LO 

CXI 

o 

CO 

o  lo 

o 

CXI 

CO 

o 

o 

LO 

o 

xf 

00 

o 

>- 

' 

<Xi  CXI 

co 

CO 

lf)LO 

LO 

LO 

X  CO 

00 

cn 

cn 

cn 

X 

o 

ro 

LO 

CT> 

CM 

o 

cn  lo 

LO 

o 

co  cn  lo 

o 

O 

X 

CO 

cn 

o 

o 

o 

o 

o 

X 

f~~ 

i —  CXI 

CO 

CO 

^  LO 

LO 

LO 

X 

X 

00 

IX 

o 

o 

O 

o 

o 

O  CO  ^ 

o 

LO 

o 

CM  OQ  Of 

cn 

O 

LO 

CM  00 

CO 

o 

cn 

CO 

X 

o 

»—  CXJ 

CO 

CO 

^  lo 

LO 

LO 

X 

X  00 

00 

cn 

cn 

o 

1 

X- 

0 

2 

LO 

00 

_ 

o 

rx  no 

cn 

o 

- —  IX  CO 

CO 

O 

LO 

r__ 

X 

CXI 

o 

o 

o 

o 

o 

1 

. —  cxi 

CXI 

co 

LO 

LO 

LO 

X 

X 

oo 

X 

0 

0 

o 

o 

O 

CO 

LO  CXJ 

00 

o 

IX 

O  LO  CXI 

o 

o 

LO 

o 

LO 

00 

CM 

LO 

o 

1 

1 —  CXJ 

CM 

CO 

lo 

LO 

LO 

X 

X 

CO 

CO 

cn 

cn 

C\J 

co 

LO 

LO 

7 

8 

cn 

o 

cxj  co  *3- 

LO 

LO 

X 

00 

cn 

o 

■ 

CXI 

CO 

LO 

f — 

1 — 

* —  r—  r— 

i — 

r - 

f— 

t — 

CXI 

CXJ 

CXI 

CXI 

CXJ 

CM 

73 


e  o  c  o  c 
o  o  o  o  c 


c  ooco-cccoooooo 

C  C  O  O  C  0  C  O  03000000 
OOOCOOOO  OOOCOOOO 
."“J  M  (M  C'i  N  CM  O'*  CN  CN  CN  CN  CN  CN  CN  CN 


o  o  o  o  o  o 


n  cn  cn  cn  rs<  cn  cn  cn  cn  cn  cn  cn  cn  cn  cn  cn 


nO 

r-» 

00 

CT\ 

j 

CN 

ro 

<r 

vC 

r- 

p»«. 

00 

r-< 

CN 

CO 

<r 

o 

sC 

O 

rH  O 

i 

i— > 

rH 

f— 4 

(—4 

, — , 

i-H 

rH 

rH 

CN 

CN 

CN 

CN 

CO 

o 

ro 

O 

O 

u* 

pH 

i-H 

c. 

X 

LO 

c 

rH 

> 

u; 

il 

li 

il 

II  II 

"C 

H 

** 

ui 

r-» 

00 

<7» 

O 

CN 

nr 

lT 

r-» 

00 

oc 

O' 

CN 

CO 

NT 

uO 

O 

1 

rH 

•—4 

rH I 

i— i 

, — i 

r—i 

rH 

i— H 

r— , 

CN 

CN 

CN 

CN 

W 

c 

U-i 

T, 

J 

s 

; — . 

o 

CQ 

H 

X 

CO 

a 

tb 

V. 

2 

< 

X 

c 

< 

£ 

lh 

X 

— 1 

CO 

,  •; 

►— ; ; 

CN 

m 

<r 

m 

r^ 

00 

o 

CN 

m 

co 

r*-» 

OC 

o 

H 

i-H 

Lt4 

! 

rH 

H 

r— ♦ 

rH 

rH 

t—4 

1— 1 

rH 

CN 

< 

< 

IQ 

U_ 

I  > 

CO 

tN 

Cl 

X 

O' 

X 

1— « 

_ 

c 

UU 

0 

CO 

< 

2 

1—1 

CO 

2 

{x. 

H 

M  w 

D-i 

o 

c 

0  0 

i— < 

i — < 

CN 

co 

<r 

\C 

r*» 

ao 

ON 

i 

CN 

CN 

ro 

r-. 

00 

O' 

O 

X 

0  < 

rJ 

i 

rH 

rH 

, — i 

rH 

rH 

rH 

rH 

rH 

O' 

C— i 

U- 

M 

w 

w 

CQ  O 

H-t 

a 

CQ 

H 

r 

o 

H 

£ 

ac 

Ch 

tH  CO 

O 

2 

S 

f— * 

O 

0 

fe 

< 

X 

0 

0 

1— » 

X 

X  X 

O 

c_- 

ari 

30 

•—i 

i-H 

UJ  H 

O 

X 

w 

On 

ag 

O' 

p 

X 

f.  §  g 

X 

x 

S. 

rH 

CN 

LO 

vO 

r--* 

00 

ON 

o 

i — i 

CN 

rH 

Hf 

UO 

vC 

< 

30 

23  0 

— 

w 

rH 

r- 

rH 

H 

4—1 

*— i 

i — 1 

w 

H 

PQ 

X 

X  DC 

J 


NODE.  LOAD  X-  Y-  Z-  X-  Y-  Z- 

NUMBER  CASE  TRANSLATION  TRANSLATION  TRANSLATION  ROTATION  ROTATION  ROTATION 


X 

x 

X 

<3- 

co 

CO 

Mt 

Mt 

vr 

mt 

X 

Mt 

MT 

mT 

MT 

mt 

X 

MT 

X 

mT 

O 

o 

o  o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

ljL 

x 

X 

X 

:  t  i 

L-J 

LU 

1 

LlJ 

1U 

LU 

LU 

LU 

LU 

X 

X 

X 

LU 

i  ,  ; 

X 

« - 

on 

CM 

o 

-X 

X 

CO 

i — 

x 

r-» 

i — 

LiO 

i — 

MT 

ON 

X 

O  X 

X 

LO 

r— 

Mt 

CC 

CM 

X  e 

1 — 

x 

X 

X 

r - 

mT 

Ox 

MT 

on 

X 

r-» 

o 

X 

X 

x 

x 

CM 

CO 

( — 

o 

c 

X 

1 — 

i — 

r— 

X 

X 

X 

X 

r-s. 

X 

X 

Mt 

CM 

X 

CO 

X 

mt 

X 

X 

X 

CM 

o 

X 

ON 

MT 

O 

X 

CM 

X 

Mf 

X 

X 

o 

LO 

CO 

1 

X 

O  1 

i 

CO 

1 

Mt 

1 

ON 

o 

1 

MT 

X 

o  / 

X 

1 

ON 

l 

ON 

1 

X 

O  I* 

X 

1 

i 

ON 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Mt 

X 

X 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

LU 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ON 

Mt 

X 

Mr 

CM 

X 

c- 

X 

X 

X 

X 

x 

X 

X 

ON 

MT 

X 

1 — 

o 

X 

X 

Mt 

X 

1 — 

CM 

X 

Mt 

1 — 

CM 

X 

X 

X 

1 — 

X 

o 

X 

X 

X 

X 

Mt 

1 — 

CNJ 

X 

1 — 

Mt 

o 

X 

o 

X 

CM 

X 

X 

Mt 

1 — 

X 

X 

1 — 

X 

ON 

o 

ON 

X 

CNJ 

X 

X 

CM 

Mt 

X 

IM* 

X 

X 

X 

ON 

CM 

X 

X 

O  1 

l 

i 

l  O  1 

1 

1 

1  O  1 

1 

1 

\  O  1 

i 

1 

l  o  o  o  o  o 

X 

CM 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

lj 

X 

t  ;  \ 

X 

X 

X 

X 

X 

!  !  , 

X 

X 

X 

X 

1 — 

X 

X 

o 

CM 

X 

X 

o 

X 

TM 

ON 

X 

X 

r— 

o 

X 

Mt 

X 

O 

Mt 

X 

X 

X 

CD 

X 

X 

X 

o 

X 

X 

f"s. 

CM 

Mt 

X 

X 

X 

X 

o 

X 

X 

X 

x 

X 

r-^ 

r— 

r— 

X 

X 

X 

(Mv 

o 

X 

t — 

ON 

O 

X 

O 

X 

X 

CD  CD  CD  CD  CD 

CM! 

CM 

o 

Mt 

Mt 

X 

X 

o 

X 

X 

X 

o 

X 

X 

X 

o 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

o 

o 

o 

o 

o 

Y 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

X 

X 

,  ,  ? 

X 

X 

X 

! - 

o 

X 

1 — 

X 

X 

X 

F - 

o 

o 

X 

Mt 

r-*. 

r>>v 

X 

Mt 

Mt 

X 

o 

X 

Mt 

cr 

X 

o 

X 

X 

Mt 

X 

i — 

i — 

X 

X 

r^ 

1 — 

X 

1 — 

X 

Mt 

o 

ON 

Mt 

o 

X 

X 

X 

X 

X 

X 

r^. 

ON 

r-v 

o 

X 

X 

r — 

X 

X 

O'. 

X 

o 

>— 

o 

X 

X 

ON 

X 

o 

X 

CO 

o 

X 

o 

ON 

ON 

Mt 

o  o  o  o  o  o 

X 

X 

X 

X 

Mt 

Mt 

X 

X 

Mt 

Mt 

Mt 

X 

Mt 

Mt 

Mt 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

X 

X 

r— 

Mt 

, — 

o 

X 

'•'•V 

, — 

, — 

ON 

CO 

CO 

ON 

o 

i — 

CD 

f— 

CO 

CO 

» — 

o 

CNJ 

CO 

o 

o 

CM 

oc 

ON 

X 

X 

X 

X 

ON 

r— 

X 

X 

CM 

ON 

CD 

X 

x 

CM 

ON 

Mt 

Mt 

X 

O' 

Mt 

X 

CNJ 

1 — 

Mt 

00 

X 

1 — 

00 

ON 

X 

CM 

X 

X 

Mt 

CM 

*— 

•— 

ON 

Mt 

•“ 

•— 

r— 

X 

•— 

•— 

X 

0  O  O  O  O  i 

i* 

1 

1 

o  » 

t 

i 

1 

O  l 

1 

1 

l’ 

O  i 

1 

i 

i  o 

X 

Mt 

X 

X 

X 

Mt 

X 

X  X; 

Mt 

X 

X  X 

mt 

X 

o 

o 

o 

o 

o 

o 

o 

o  o 

o 

o 

o 

o 

o 

o 

o 

X 

X 

X 

X 

X 

X 

X 

X  X 

X 

X 

X  x 

X 

X 

X 

o 

O 

ON 

CD 

o 

X 

1 — 

x  o 

1 — 

X 

^  ON 

CO 

X 

mT 

X 

i — 

o 

t— 

CO 

X 

CO 

<—  00 

o 

X 

1 —  X 

CD 

CO 

X 

X 

ON 

CNJ 

X 

X 

X 

00 

Mt  O 

00 

Mt 

Mt 

Mt 

X 

CO 

X 

X 

ON 

Mt 

Mt 

o 

X 

X 

r—  ON 

Mt 

o 

Mt 

X 

o 

Mt 

X 

Mt 

CM 

f - 

X 

X 

X 

1 — 

x  • — 

X 

1 — 

X 

X 

1 — 

X 

• — 

O  i 

O  1 

’  O  i# 

O  i 

o  o  o  o  o 

in  n  c\j  i —  ocn'XiN'LDLr)^roc\jr-oa'OON^Ln^ro(\j 

CM  CM  CNJ  CNJ  CNJ  CM  i —  I —  r—  i —  * —  ■ —  i —  i —  r—  r — 


1 558L  t-02  . 1698E+02  .2221EM1  . 3195E  +  Q1  .7731E+00  -.4641E-01 


O  CO 

- —  LO 

• —  v? 

o 

o  o 

O  O 

o  o 

O  O 

!  ,  ■ 

LU  LU 

UJ  LU 

LU  U 

LU  s_U 

cm 

uD  CO 

« —  i — - 

LO  oo 

LO 

CM  oo 

CM  OO 

CM  CD 

CO  P". 

' 

cp 

O  LO 

cp  o 

MT  vO 

CD  L D 

j 

VO 

1 

1  i 

LO  CM 

1  1 

CM  CD 

CM  CD 

co 

< —  lo 

o  mt 

O  MT 

O  CD 

\ 

o 

o  o 

o  o 

O  O 

O  O 

i 

i 

i  i 

■+■  1 

4-  1 

4-  l 

i 

UJ 

UJ  UJ 

LU  LU 

LU  Lj 

LU  LU 

4 

O 

O  i— 

i —  p^. 

CD  LO 

LO  00 

lo  i — 

i —  f\ 

cm  r^. 

CM 

lO 

00  LO 

LO  O 

00  LO 

O  00 

1 

CM  co 

1 —  cm 

CM  CD 

CP  r— 

iD 

r—  MT 

i -  LO 

< —  OO 

. —  CD 

O 

o  o 

O  O 

o  o 

O  O 

i 

4-  1 

+  1 

4-  1 

+  1 

LU 

LU  LU 

LU  LU 

LU  LU 

LU  LU 

CM 

CD 

O  O 

LO  LO 

CM  LO 

LO  r— 

r^. 

CP  LO 

CM  r- 

r- 

* —  LO 

CO  lO 

mT  CP 

LO  CP 

i — 

r—  LO 

CM  i — 

CD  <— 

CD  > — 

• -  <P 

) — 

• —  mT 

O  LO 

O 

o  o 

o  o 

o  o 

o  o 

1 

+  1 

4-  1 

+  1 

4-  » 

i_U 

LU  LU 

LU  LU 

LU  Uj 

LU  LU 

t — > 

o  o 

r\ 

CP  LO 

CP  O 

i — 

oc  pt 

oo  O 

CM  LO 

P*^  mT 

r— 

CM  i— 

LO  1-"-. 

1—  o 

^  CM 

< - 

mt  cm 

CD  CM 

CM  « — 

LO  CD 

CM  lO 

, — 

LO 

CM 

LO 

CM 

LO 

CM 

LO 

O  O 

o 

o 

o 

o 

O 

o 

O 

o 

■f  1 

-*■ 

1 

+ 

1 

+ 

1 

4- 

l 

U  LU 

•U 

LU 

LU 

LU 

uJ 

LU 

UJ 

■u> 

00  o 

vr 

i — 

CP 

CP 

r— 

LO 

1 — 

y 

CP  CD 

CO 

00 

CP 

co 

o 

CO 

O 

LO  oo 

■M" 

LO 

o 

CP 

LO 

LO 

rr. 

LO 

1 — 

CD 

CP 

r“ 

CM 

' — 

1 — 

JO 

CD 

CO 

CM 

lO  LO 

0 

O 

O 

O 

o 

"xT 

O  Mp 

LO 

0 

O 

O 

O 

_ 1 

LO 

, — 

0 

O 

O 

O 

CM  lO 

i — 

LO 

CM 

LO 

CM 

LO 

CM 

LO 

o  o 

O 

o 

o 

o 

O 

o 

O 

o 

l 

4- 

1 

4- 

1 

-I- 

1 

4 

1 

LU 

.  U  LU 

LU 

!  i  : 

LU 

LU 

LU 

LU 

LU 

LU 

00  CD 

< — 

CD 

CP 

O 

LO 

CP 

o 

5: 

II 

II 

II  II 

II 

11 

II 

II 

II 

LO  CO 

00 

LO 

o 

CP 

o 

o 

LO  CD 

CM 

i-— 

CD 

LO 

co 

CP 

LO 

♦— * 

CM 

r—  MP 

r— 

i — 

CD 

i — 

mT 

« — 

lO 

o 

i — 

CO 

zz 

i — 

* 

o 

O0 

r— 

k-m 

2 

h - < 

CO 

1 — 

O 

O0 

OO 

2 

2 

—< 

U- 1 

>- 

O 

r—  r— 

r - 

i — 

r— 

r— 

( — 

« — 

i — 

r - 

o 

2: 

h- 

O0 

_U 

►— 1 

o 

n 

ti 

I: 

OC 

c 

2 

>- 

< 

h- 

_ 1 

o 

5: 

O 

_ 1 

2 

C 

> — ) 

i — 

LU 

QC 

*— • 

«=t 

<C 

OC 

o 

h- 

0 

I — 

2 

Q_ 

LU 

1 — 

OO 

Lu 

CJ 

<3: 

^r 

UJ 

o 

h- 

>- 

ZD 

OO 

h~ 

C 

ZD 

h- 

:o 

• — 1 

ZD 

QC 

2: 

Q_ 

LU 

ZD  0O 

O0 

or 

LU 

OC 

2 

h— 

O 

LU 

2C 

2C 

O-  OO 

*— % 

1 — 

OO 

1 — 

» — < 

_ I 

ZD 

L* 

1 — < 

► — i 

Lu 

2:  lu 

OO 

X 

2: 

c_> 

CM  CM 

CD 

CD 

LO 

LO 

LO 

LO 

_ 1 

h- 

O 

Ll. 

*— «  2: 

LU 

O 

LU 

CL 

r—  r— 

< — 

i — 

i — 

r— 

< — 

i — 

i — 

f— 

o 

O 

o 

h- 

\— 

i — • 

U_ 

_ 1 

Q_ 

Cl 

LU 

OO 

LU 

LU 

h~ 

CD  Lu 

c 

LU 

OO 

00 

| — 

LO 

OC 

» — i 

O0 

<C  *— « 

2 

ZD 

LU 

00 

2 

LU 

_ 1 

O 

OH< 

J 

DC 

LU 

1 

o 

CO 

oo 

Cl. 

H- 

_J  oO 

<C 

00 

>- 

CO 

> — < 

<C 

oo 

CJ 

Q 

2 

CD 

h— 

— 1 

oo 

_ 1 

LU 

_J  _i 

1— 1 

2: 

LU 

O 

l 

♦— » 

c 

O- 

UJ 

<x. 

<C 

h- 

LU 

CO 

CL 

Cl 

ZD 

UD 

cc 

Q 

LU 

O  1— 

O 

QC 

OO 

UJ 

1— 

cr 

* — < 

h- 

QC 

O 

_J  o  o 

1— 

t— 1 

O 

LU 

h- 

LU 

o 

oo 

z 

LU 

2C 

00 

LU 

Lu 

CL 

OO 

<c  LU 


iD 


OO 


o 


79 


STEP-BY-STEP  INTEGRATION 


SAP  4A  ’  EST  CASt  CURVFl)  PL 


ro 

r- 

/M 

LO 

ro 

r>* 

L.0 

ro 

fx. 

lO 

ro 

r>. 

r*^ 

lO 

ro 

r^. 

LO 

ro 

XT 

lo 

cn 

LO 

xr 

lo 

O' 

LO 

*3" 

LO 

cr^ 

lO 

LO 

CX 

LO 

ro 

LO 

o> 

+ 

+ 

(XI 

xr 

lO 

lo 

CXI 

LO 

LO 

CXI 

LO 

LO 

CM 

xf 

LO 

LO 

CM 

LO 

LO 

lO 

lO 

o 

o 

o 

O 

o 

lO 

lO 

vO 

LO 

CM 

LO 

uO 

LO 

LO 

lO 

LO 

lO 

cO 

lO 

LO 

LO 

CXI 

CXl 

CXi 

CM 

CXi 

uo 

LO 

lO 

LO 

LO 

rv 

r— 

i-"*. 

r». 

O 

o 

o 

o 

o 

o 

o 

o 

O 

r~- - 

CXI 

cm 

CXi 

CXI 

CM 

ro 

ro 

ro 

ro 

ro 

lO 

LO 

LO 

LO 

LO 

o 

lO 

lO 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

ro 

o 

CM 

LO 

r^. 

O 

CM 

lO 

O 

CM 

LO 

rM 

O 

CM 

LO 

O 

CM 

lO 

r— 

O 

4* 

ro 

4- 

O'. 

o  1— 

CM 

ro 

.-O 

O  r— 

CM 

ro 

lO 

o  *— 

CM 

ro 

cO 

O  1— 

CM 

ro 

lO 

o  r— 

CM 

ro 

lO 

rn 

no 

cn 

ro 

LO 

ro 

CM 

o 

— - 

, _ , 

, _ . 

, — 

r- 

i— 

ro 

xr 

REFERENCES 


K.C.  Bathe,  E.L.  Wilson,  and  F.E.  Peterson,  SAP  IV  -  A 
Structural  Analysis  Program  for  Static  and  Dynamic  ^Response 
of  Linear  Systems,  University  of  California,  Berkeley,  CA. 
EERC  73-11,  April  1974. 

M.A.  Pa lie,  A  Quadrilateral  Plane  Stress  Finite  Element  with 
Bending/Extensional  Coupling^,  Thesis  University  of  Lowell, 
December  1976. 


R.M .  Jones,  Mechanics  of  Composite  Materials,  McGraw-Hill 
Book  Company*  New  York,  197  5  .  .  — 


